Disentangling multipole resonances through a full x-ray polarization analysis. 
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Complete polarization analysis applied to resonant x-ray scattering at the Cr K-edge in K2Cr04 
shows that incident linearly polarized x-rays can be converted into circularly polarized x-rays by 
diffraction at the Cr pre-edge [E = 5994 eV). The physical mechanism behind this phenomenon is a 
subtle interference effect between purely dipole (El-El) and purely quadrupole (E2-E2) transitions, 
leading to a phase shift between the respective scattering amplitudes. This effect may be exploited 
to disentangle two close-lying resonances that appear as a single peak in a conventional energy scan, 
in this way allowing to single out and identify the different multipole order parameters involved. 

PACS numbers: 78.70.Ck, 78.20.Bh, 78.20.Ek 



I. INTRODUCTION 

In the last 10 years resonant x-ray scattering (RXS) 
has developed into powerful technique to obtain direct in- 
formation about charge, magnetic, and orbital degrees of 
frccdoroii^ii^ii^i^i^. It combines the high sensitivity of x-ray 
diffraction to long-range ordered structures with that of 
x-ray absorption spectroscopy to local electronic config- 
urations. In particular, the development of third genera- 
tion synchrotron radiation sources has made possible the 
detection of small effects in electronic distribution, due 
to magneto-electric anisotropy^ or to local chiralit3*2i>^, 
that can be related to the interference between dipole 
(El) and quadrupole (E2) resonances. These pioneer- 
ing studies paved the way to a new interpretation of 
RXS experiments in terms of electromagnetic multipoles 
of higher order than dipole, and led to the detection 
of phase transitions characterized by order parameters 
(OP) of exotic symmetryiSiiiii^. Several theories have 
been developed, based on these higher-order OP, to ex- 
plain "anomalous" phase transitions. For example, in 
Np02 a proper interpretation of the magnetic ground- 
state requires a primary OP at least octupolar order—, 
whereas in high-temperature cuprate superconductors, 
the pseudogap phase has been interpreted in terms of 
parity and time-reversal odd toroidal multipolesi^ii^. In 
several cases, though, the assignment of the multipolar 
origin to a RXS signal is not clear— li^iii. The char- 
acteristic variation of the intensity and polarization as 
the sample is rotated about the scattering vector dur- 
ing an azimuthal scan may allow a clearer identification 
of the order of the multipole. However, this technique is 
plagued by many experimental difficulties either from the 
sample, e.g. when the crystal presents twinning and mo- 
saicity, or due to restrictive sample environments, such 
as cryomagnets. Moreover, it is very difficult to iden- 
tify and resolve two resonances of different multipolar 



origin when they are separated by less than ~ 1 eV, a 
situation which frequently occurs at the metal pre K- 
edge region of transition metal oxides, where El and 
E2 transitions can have similar magnitudes, or at L2.3 
edges of rare-earth compounds. For example, at the Fe 
pre K-edge in a-Fe203, evidence of a E2-E2 transition 
was found already in 1993^^. This transition was inter- 
preted as the signature of an hexadecapolar electronic 
orderingi^. Later, however, it was shown by symme- 
try arguments that also an axial-toroidal-quadrupole OP, 
parity-breaking, was hidden within the same resonanceii, 
but to date no conclusive experimental evidence for this 
interpretation has been provided. Analogously, at the 
Ce L2,3-edges in CeFe2 the different electronic OP re- 
lated to 4/ and 5d states are entangled and cannot be 
examined individually by conventional RXS^. A further 
example is provided by K2Cr04. Its space group symme- 
try [Prima, No. 62) allows several excitation channels at 
the (IkO) Bragg forbidden reflections^, due to the pres- 
ence of multiple moments of different symmetries (elec- 
tric quadrupole, octupole, hexadecapole) in the same en- 
ergy range, as described in more details below. 



The aim of the present article is to address the previous 
problems by extending the well-known techniques of op- 
tical polarimetry from the visible to the x-ray regime, as 
developed at the beamline ID202^ at the ESRF, Greno- 
ble, France. By using a diamond x-ray phase plate to ro- 
tate the incident linear polarisation in combination with 
a linear polarisation analyser, we can resolve resonances 
determined by multipoles of different order that are very 
close in energy, playing on their relative phase shifts. The 
idea can be explained through a simplified model, by con- 
sidering two externally driven, damped harmonic oscilla- 
tors of unitary amplitude, but with resonant frequencies 
differing by 2^. The scattering amplitude of such oscilla- 
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tors is given by 



uj±C~iT 



(1) 



Here r is the inverse damping time, and uj the frequency 
of the external excitation. We also suppose that the two 
resonances scatter in different polarization channels. Us- 
ing the Jones matrix formalism^^ the polarisation of the 
scattered beam may then be written as e' = Ge, where e 
(e') arc the incident (scattered) polarisation vectors, and 

9- 9+ 
9+ 9- 

the photon energy. Experimentally, the scattered beam 
polarisation is best described in terms of the Poincare- 
Stokes parameters: 
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with Pq = (|ejjp-|- |e^p) the total intensity, and e'* the 
complex conjugate of e'. P{ and P2 describe the linear 
polarization states, whereas P3 indicates the degree of the 
circular polarization. The Poincare-Stokes parameters, 
^0,1.2,3 of the incident beam arc obtained by substituting 
e for e' . 

For example, for an incoming tt polarised beam, e = 



wc obtain 



, which in turn yields: 
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Therefore we expect that in the intermediate region be- 
tween the two resonances, the outgoing beam is circularly 
polarized, depending on the relative dcphasing of the two 
resonances, as shown in Fig. [1] for the case C = T. 

Below, we describe experimental data which we then 
compare to quantitative ab-initio calculations carried out 
using the FDMNES code^. We demonstrate that a 100% 
linear- to circular-polarization conversion at the prc-edge 
region of the Cr K-cdge in K2Cr04 is induced by the in- 
terference of the dispersive and absorptive parts of two 
different multipolcs probed by purely dipole (El-El) and 
purely quadrupole (E2-E2) resonances. Thus, the scat- 
tered beam originates from two different excitation chan- 
nels, each scattering the beam with a different phase. 

Their relative amplitudes, at a given energy, arc gov- 
erned by the probed multipoles, the relative orientation 



FIG. 1: (Color) - Top: Typical behavior of real {g' , continu- 
ous) and imaginary {g" , dashed) parts for the two resonators 
of Eq. [1] {g+ red; g_ blue), with ( = T. Middle: Corre- 
sponding phase shifts. Bottom: Poincare-Stokes parameters 
as calculated from Eqs. |3] 



of the incident and scattered electric field vectors to the 
probed multipoles and the reciprocal lattice point under 
study. Their interference can therefore be tuned simply 
by varying the incident polarization by means of a phase 
plate. 



II. EXPERIMENTAL SETUP 

Experiments were carried out at ID20, ESRF. The ex- 
perimental setup is outlined in Fig. [51 A single crystal of 
K2Cr04 was mounted on the six-circle horizontal diffrac- 
tomcter, and a cryostat was used to stabilize the temper- 
ature at 300 K. Sample rocking curves (6'-scans) of the 
crystal resulted in widths smaller than 0.01° indicating 
a high sample quality. 

A diamond phase plate of thickness 700 |^m with a 
[110] surface was inserted into the incident beam, within 
its own goniometer, and the (111) Bragg reflection in 
symmetric Laue geometry was selected to modify the po- 
larization of the x-ray beam incident on the diffractome- 
ter. The phase plate was operated in either quarter- wave 
or half-wave plate mode. With the former we generated 
left- or right-circular polarization, P3 « ±1, whereas with 
the latter we rotated the linear incident polarization into 
an arbitrary plan o^^i^^ , described by Pi « cos(277) and 
P2 ~ sin(27/). Here ry is the angle between the inci- 
dent beam electric field vector and the vertical axis (see 
Fig. [2]), i.e. 77 = when the polarization is perpendicular 
to the horizontal scattering plane (cr polarization). 

The sample was mounted with the [010] and [100] di- 
rections defining the horizontal scattering plane. Fig- 
ure [3] shows the fluorescence yield and the energy depen- 
dence of the glide-plane forbidden (130) reflection, col- 
lected at an azimuthal angle of 5* = —0.78(3)° degrees 
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FIG. 2: (Color online) Experimental setup with phase plate 
in half-wave mode; x-ray directions indicated by blue arrows, 
polarizations by red ones. Synchrotron light arrives from the 
left, horizontally polarized (tt). is the rotation angle of the 
incident polarization, rj' is the rotation angle of the polariza- 
tion analyser crystal; the zero positions of the two angles, cor- 
responding to a and a' polarisations respectively, are repre- 
sented by dashed lines. The continuous line is the the rocking 
axis of the polarisation analysis crystal 9pA- The polarisation 
analyser stage is shown in the configuration corresponding to 
the maximum detected intensity. 
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FIG. 4; (Color online) Variation of the integrated intensity 
of the (130) forbidden Bragg peak as a function of the polar- 
ization analyzer angle rj' . The dashed line represents a fit to 
these data with Eq. [3] (see text for details). The insert is an 
example of a rocking scan of the analyser: the integrated in- 
stensity obtained from the fit (dashed line: lorentian squared 
and linear background) represents experimental point in the 
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intensities were then fitted to the equation 
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FIG. 3: (Color online) Experimental data on fluorescence 
yield (black continuous line) and energy scan (red symbols) 
for the (130) reflection in tt — > a' polarization configuration. 
Arrows indicate the energy values where Stokes' parameters 
were measured. 



with respeet to the reciprocal lattice direction [010]. The 
photon energy was then tuned to the pre-edge region of 
the Cr K-edge (5994 eV). An x-ray polarization analyzer 
was placed in the scattered beam. It exploited the (220) 
Bragg reflection of a LiF crystal, scattering close to Brew- 
ster's angle of 45°. The polarization analyzer setup was 
rotated around the scattered beam by an angle, 77', and 
at each point the integrated intensity was determined by 
rocking the analyzer's theta axis (^pa). An example is 
shown in the inset of Fig. 01 The resulting integrated 



P' 

I = ^[l + P[ cos 2ri' + sin 2r]'] , (3) 

to obtain the Poincare-Stokes parameters, P^ and Pj- An 
example is shown in Fig. [5] for 77 = -90°. The degree of 
circular polarization, P3, can not be measured directly 
in this setup. However, an upper limit is inferred from 
+ + < 1 (tlie equality holding for a completely 
polarized beam). 

We systematically measured P{ and P2 of the beam 
scattered at the (130) reciprocal lattice point as fimc- 
tion of jy. Figure O shows both the experimental data 
(symbols) and the theoretical calculation (dashed lines), 
described below. The measured degree of linear polariza- 
tion of the scattered beam, P{^ + P!^ ^ strongly deviates 
from 100% in the range —10° ^ ij ^ 50°, indicating that 
a large component of the scattered beam is either cir- 
cularly polarized or depolarized. To ascertain which of 
these two processes is realised we reconfigured the dia- 
mond phase plate to produce circularly polarized x-rays. 
Figure [S] shows the measured P{ and Pj for both lin- 
early and circularly polarized x-rays. The presence of 
linearly scattered x-rays for the circular incident case is 
consistent with the assumption that the region for which 
Pl^ + P2 strongly deviates from 1 corresponds to an in- 
creased circularly polarized contribution. Furthermore, 
the calculations for P^ , detailed in Section III, involving 
only completely polarized contributions to the scattered 
beam, agree well with its upper limit inferred from the 
data (Pg^ = 1 — P!^ — -P2^)j indicating that the signal is 
essentially circularly polarised. 
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FIG. 5: (Color online) Calculated (dashed lines) and mea- 
sured (symbols) Stokes' parameters for (130) reflection versus 
the polarisation angle of the linear incoming light (see Fig. [2| , 
E = 5.994 keV. See text for details. 
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FIG. 6: (Color online) Stokes' parameters for (130) reflection 
with (top) linearly (77=10°) and (bottom) circularly polarized 
incident x-rays, collected at iJ = 5.994 keV. Dashed lines are 
flt with Eq. [3] 



III. THEORETICAL DISCUSSION 

In RXS the global process of photon absorption, vir- 
tual photoelectron excitation and photon re-emission, is 
coherent throughout the crystal, thus giving rise to the 
usual Bragg diffraction condition J2j g'*^ ' ifoj + /j + 



if'j)- Here Rj stands for the position of the scatter- 
ing center j, Q is the diffraction vector and foj is the 
Thomson factor, /j and related by Kramers-Kronig 
transform, are given, at resonance, by the expression^: 



{Mj)\0'*\^n){i^n\0\Mj)) 
UJ - {uJn- tOg) - 1 



(4) 

where uj is the photon energy, cOg the ground state en- 
ergy, LUn and r„ are the energy and inverse lifetime of 
the excited states, tpg{j) is the core ground state centered 
around the j*'' atom and tpn the photo-excited state, e 
and e' are the polarizations of the incoming and outgoing 
photons and q and q' their corresponding wave vectors. 
The sum is extended over all the excited states of the sys- 
tem. The transition operator ^ = e^^ ■ r (l — ^ef ^ ■ r) 
is written as a multipolar expansion of the photon field 
up to electric dipole (El) and quadrupole (E2) terms; r is 
the electron position relative to the resonating ion, •* is 
the polarization of the incoming (outgoing) photon and 
' its corresponding wave vector. 

By taking into account the space group symmetry of 
K2Cr04 {Pnma, No. 62), the four equivalent Cr sites at 
Wyckoff 4c positions (with local mirror-plane rhy) can 
be related one another by the following symmetry oper- 
ations: /a = //i, /4 = C'2a;./i^and /2 = 1/4. / is the 
space-inversion operator and C'2x is the 7r-rotation op- 
erator around x-axis. The Thomson scatering foj does 
not contribute at (IkO) type reflections, for any k, due 
to the glide plane extinction rule. Therefore, the struc- 
ture factor at Cr K-edgc for the (IfcO) chosen reflections, 
when summed over all equivalent sites, becomes, for El- 
E2 scattering: 



k 

S = 2i sin[27r(a; + -)](! - m^)/! 
while for El-El and E2-E2 scattering it is given by: 



S = 2cos[27r(x-h -)](1 



(5) 



(6) 



Here x ~ 0.23 is the fractional coordinate of Cr atoms 
and rhz is a glide-plane orthogonal to the z-axis. In de- 
riving Eqs. [5]and[6l we have used the identity /i — rhyfi. 
It is interesting to note the different behavior of the two 
terms for k even or odd. For example, when fc = 4, the 
E1-E2 scattering is proportional to sin(27ra;) ~ 0.99 and 
it dominates the other terms, proportional to cos(27ra;) ~ 
0.12. Indeed we found the presence of a very intense 
pre-edge feature from E1-E2 channel at the (140) re- 
flection that is related to the electric octupolc moment, 
as predicted in Ref. [l^ and verified numerically by our 
ab-initio calculations. However, as described above, the 
presence of a single scattering channel, as in the case of 
the (140) reflection, can not lead to a circularly polarized 
diffracted beam. This can be demonstrated by a sym- 
metry argument: if only one scatterer is present, which 
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is by hypotesis non-magnetic, and the incident light is 
Unearly polarized, then the initial state is time-reversal 
even. Therefore, as matter-radiation interaction docs not 
break time-reversal, it follows that the final state must 
also be time- reversal even, i.e., radiation can not be circu- 
larly polarized, which would break time-reversal. This is 
no more true when two scattcrers arc present, due to the 
extra degree of freedom represented by the time (phase) 
delay between the two scattering processes. Indeed, this 
was a posteriori verified by our numerical simulation, 
which confirmed that no outgoing circular polarization 
is present at the (140) reflection. 

The case of the (130) reflection is very different. The 
role of sin(27rx) and cos(27ra;) in Eqs. [5] and [6] switches 
in such a way that the two diffraction channels El-El 
and E2-E2 become predominant. Further analysis of 
the structure factor— reveals that two resonances arc al- 
lowed, one for each channel, corresponding to an elec- 
tric quadrupole ordering for the El-El scattering and 
an electric hexadecapolei^ for the E2-E2 scattering. Fi- 
nally, multiple-scattering calculations with the FDMNES 
program confirm that the two resonances overlap in the 
pre-edge region, though slightly shifted in energy of ~ 1 
eV. These are the conditions to be met to get the interfer- 
ence of the two channels. In order to describe the effect 
quantitatively from a theoretical point of view we used 
the ab-initio code FDMNES, in the multiple-scattering 
mode, to calculate P3 directly from Eq. |4]for the K2Cr04 
structuro2i. We employed a cluster of 43 atoms, cor- 
responding to a radius of 5.5 A around the resonating 
Cr-atom. Notice that in this most general case we find 
P^iio) (X (/^i(o.)/i^2(^)7/B2H/Bi(^)) where /' and 
/" are the usual dispersive and absorptive terms (see 
Eq. 3]) for El and E2 channels. Therefore at the photon 
energy w, Pg is determined by the interference of the ab- 
sorptive part /" of one channel with the dispersive part 
/' of the other and, again, in the presence of just one 
channel, P3 = 0. 

The numerical simulations shown in Fig. O (dashed 
lines) confirm that there is an azimuthal region where 
the incoming linear polarization is fully converted into 
an outgoing circular polarization and that the effect is 
determined by the interference of the El-El and E2-E2 
channels. 

In this respect, it is interesting to note that the ori- 
gin of this effect is profoundly different from those de- 
termined by a chiral (magnetic) structure (see, e.g., Ref. 
[29ll30h . as clearly seen by the fact that all the tensors 



involved are non-magnetic and parity-even. 

Finally, we verified experimentally that at the Cr K- 
edge {E = 6010 and 6018 eV), where only one term in the 
El-El channel is present, no circular polarization was ob- 
served for all incident angles, i.e., + P!^ — 1. Calcu- 
lations performed using FDMNES confirmed this result. 
IV. CONCLUSIONS 

Polarization analysis of RXS experiments has devel- 
oped greatly in the last few years, helping to understand 
several characteristics of order parameters in transition 
metal oxides, rare-earth based compounds, and actinidcs. 
Up to now, however, the full investigation of Stokes' pa- 
rameters was not applied most likely because linear po- 
larization analysis, where only the P[ parameter was de- 
termined by measuring the a a' and a ^ it' chan- 
nels, was considered sufficient. While this may be true 
when just one excitation channel is involved (as at the 
(140) reflection in the present case), several dephasing 
phenomena may appear when two different multipolc ex- 
citations close in energy are involved in the transition. 
As we have seen, these phenomena may lead to a situa- 
tion where incoming linear polarization is scattered to a 
circular polarization due to an interference between two 
multipoles, at the same time allowing for a very sensi- 
tive determination of the presence of the second multi- 
pole. We believe that the use of phase plates and of a 
complete polarization analysis, is the key to disentan- 
gle multi-resonance structures in those situations where 
an usual energy scan, like the one shown in Fig. [3l is 
not sufficient to this aim. Effect of d-band filling on de- 
tails of the electronic structure will be investigated by 
the method presented here in the series of isostructural 
compounds K2Cr04 K2Mn04 ^ K2Fe04. 
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